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Parity-odd domains, corresponding to non-trivial topological solutions of the QCD vacuum, might 
be created during rclativistic heavy-ion collisions. These domains arc predicted to lead to charge 
separation of quarks along the orbital momentum of the system created in non-central collisions. 
To study this effect, we investigate a three particle mixed harmonics azimuthal correlator which is 
a T'-even observable, but directly sensitive to the charge separation effect. We report measurements 
of this observable using the STAR detector in Au+Au and Cu+Cu collisions at ^/snn=2QQ and 
62 GeV. The results arc presented as a function of collision ccntrality, particle separation in rapidity, 
and particle transverse momentum. A signal consistent with several of the theoretical expectations is 
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detected in all four data sets. We compare our results to the predictions of existing event generators, 
and discuss in detail possible contributions from other effects that are not related to parity violation. 

PACS numbers: ll.30.Er, ll.30.Qc, 25.75.Ld, 25.75.Nq 



I. INTRODUCTION. 

Quantum Chromodynamics (QCD) is widely accepted 
as the theory of the strong interaction. The perturba- 
tive regime, applying to processes with large momentum 
transfer, is theoretically calculable and has been exten- 
sively tested [1]. On the other hand, the regime in which 
quarks and gluons interact with modest momenta and 
with an effective coupling constant that is too large for 
perturbation theory to apply, cannot be reliably calcu- 
lated by analytic methods. Lattice gauge theory is one 
first principle approach which can be used. It predicts 
the existence of a new state of strongly-interacting mat- 
ter at high energy density. This state has now been ob- 
served in high-energy heavy-ion collisions at the Rela- 
tivistic Heavy Ion Collider (RHIC) at Brookhaven Na- 
tional Laboratory P]. 

Many interesting features of this new state of mat- 
ter produced in these collisions have been observed. In 
this paper we focus on a new phenomenon which we re- 
fer to as local strong parity {V) violation. It is well 
known that the strong interaction conserves parity — 
meaning that strong interactions do not lead to reac- 
tions which produce a finite expectation value for any 
parity odd (changing sign under parity transformation) 
observable. The best evidence for this comes from ex- 
periments which set limits on the electric dipole moment 
of the neutron [3l |4] . These experiments show that the 
parameter 9 whose nonzero value would describe parity 
violation in QCD must have magnitude less than lO^^'^. 
This limit effectively makes direct global parity violation 
unobservable in heavy-ion reactions. Our measurement 
of a similar T'-odd observable is consistent with zero at 
the level of the experimental precision of 10~^ (see sec- 
tion 



tions with non-zero topological charge have a finite ex- 
pectation value for {Echromo ■ Bchromo), whcrC Echromo 
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The concept of local parity violation at high tempera- 
tures or in high-energy heavy-ion collisions was discussed 
by Lee and Wick [51 [S] and Morley and Schmidt Ji and 
elaborated by Kharzeev et al. [5]. In a dense highly 
excited state, gluon fields can produce configurations, 
local in space and time, which cause V, time rever- 
sal T, and, via the CPT theorem, CV violating effects. 
These field configurations form in different ways in dif- 
ferent events and averaged over many events they would 
not yield a finite expectation value for a T'-odd observ- 
able. Each space-time region, occupied by such a config- 
uration, is spontaneously produced with a random sign 
of T'-violation, which in the theory is determined by 
the gluonic field topological charge^. Field configura- 



^ The topological charge distinguishes gluonic field configurations 



and Bchromo are the chromoelectric and chromomagnetic 
fields, and the average is taken over the region occupied 
by the configuration. Since the space time symmetries 
of chromodynamic fields are the same as those of elec- 
tromagnetic fields, with electric field being a vector and 
magnetic field being a pseudovector, this region is not in- 
variant under V (and T) transformations. Quark interac- 
tions with such topological gluonic configurations change 
the quark chirality leading to asymmetry in the num- 
ber of left and right handed quarks: Nj^ — Nj^ = 2nfQ, 
where nf is the number of light quark flavors and Q is 
the topological charge of the gluonic configuration. Thus, 
the gluonic field configurations with non-zero topological 
charge induce the local T'-violating effects. Different as- 
pects of an experimental detection of this phenomenon 
were discussed in [HI [T^] . 

In non-central collisions such a domain can manifest 
itself via preferential same charge particle emission along 
the system angular momentum [T31[T3] (see Fig. 1). Op- 
posite charge quarks would tend to be emitted in oppo- 
site directions relative to the system angular momentum. 
This asymmetry in the emission of quarks would be re- 
fiected in, for example, an analogous asymmetry between 
positive and negative pion emission directions. This phe- 
nomenon is driven by the large (electro-) magnetic field 
produced in non-central heavy-ion collisions [131 [TSl US] . 
Peak magnetic field strengths can reach levels of the 
order of 10^^ T. The combined effect of this magnetic 
field (which tends to align the magnetic moments of the 
quarks with the field) and the difference in the number 
of quarks with positive and negative chiralities (which 
is induced by their presence in a "7^-violating bubble") 
gives rise to the "Chiral Magnetic Effect" . 

The same phenomenon can also be described in terms 
of induction of electric field by the (quasi) static mag- 
netic field, which occurs in the presence of these topo- 
logically non-trivial vacuum solutions il6^. The induced 
electric field is parallel to the magnetic field, and leads to 
the charge separation in that direction. Thus the charge 
separation can be viewed as a non-zero electric dipole 
moment of the system (see Fig. [T]) . 

Depending on the sign of the domain's topological 
charge, positively charged particles will be preferentially 
emitted either along, or in the direction opposite to, the 



that can not be continuously transformed one into another. In 
general it is not expected to be "neutralized" and in a given event 
the net topological charge can take non-zero values. For a review 
of topological effects in QCD, see [9l I10| 



4 



LorS A 




FIG. 1: (Color online) Schematic view of the charge separa- 
tion along the system orbital momentum. The orientation of 
the charge separation fluctuates in accord with the sign of the 
topological charge. The direction of the orbital momentum L, 
and that of the magnetic field B, is indicated by an arrow. 



system orbital angular momentum, with negative par- 
ticles flowing oppositely to the positive particles. The 
magnetic field and the angular momentum are normal to 
the plane containing the trajectories of the two collid- 
ing ions. This plane, called the reaction plane, can be 
found experimentally in each collision by observation of 
the azimuthal distribution of produced particles in that 
event. 




FIG. 2: (Color online) Schematic view of the transverse plane 
indicating the orientation of the reaction plane and particle 
azimuths relative to that plane. The colliding nuclei are trav- 
eling into and out of the figure. 

When two heavy ions collide with a finite impact pa- 
rameter, the probability for particles to be emitted in 
a given azimuthal direction is often described with a 
Fourier decomposition |17j : 



where A(f) = {4> — ^rp) is the particle azimuthal direction 
relative to the reaction plane, as shown in Fig. 2. vi and 
V2 are coefficients accounting for the so-called directed 
and elliptic fiow, respectively, and a indicates the parti- 
cle type. They depend on the impact parameter of the 
colliding nuclei, the particle type (tt, K , p, ...), transverse 
momentum (pt), and pseudorapidity (77) of the produced 
particles. For collisions of identical nuclei, symmetry re- 
quires vi to be an odd function of rapidity and V2 to be 
an even function of rapidity. Measurements (for a re- 
view and references, see [TB]) have found that, at RHIC, 
vi is quite small at mid-rapidity; typically, < 0.005 
for — 1<77<-|-1. In contrast, V2 is found to be siz- 
able and positive. In Au-|-Au collisions at y/sMN= 200 
GeV, for unidentified charged hadrons, V2 reaches 0.25 
for pt ~ 3 GeV/c, and 0.06 when integrated over all pt- 

Phenomenologically, the charge separation due to a do- 
main with a given sign of the topological charge can be 
described by adding P-odd sine terms to the Fourier de- 
composition Eq. [1] [T^ : 

" cx 1 -I- 2ui_ct cos(Z\0) -I- 2 i;2.Q cos(2Z\0) + ... 



2ai^Q, sin(Z\(/)) -f 2 02,0 sin(2Z\0) 



(2) 



cx 



1 -I- 2wi,Q cos{Acf)) + 2 V2,a cos{2A(t)) + ... ,(1) 



where the a parameters describe the T'-violating effect. 
Equation [2] describes the azimuthal distribution of par- 
ticles of a given transverse momentum and rapidity and, 
like the flow coefficients, a coefficients depend on trans- 
verse momentum and rapidity of the particles. In addi- 
tion, they depend also on the rapidity (position) of the 
domain. One expects that only particles close in rapidity 
to the domain position are affected. According to the 
theory, the signs of a coefficients vary following the fluc- 
tuations in the domain's topological charge. If the par- 
ticle distributions are averaged over many events, then 
these coefficients will vanish because the distributions 
are averaged over several domains with different signs 
of the topological charge. However, the effect of these 
domains on charged particle correlations will not van- 
ish in this average, as discussed below. In this analysis 
we consider only the first harmonic coefficient ai, which 
is expected to account for most of the effect although 
higher harmonics determine the exact shape of the dis- 
tribution. For brevity we will omit the harmonic number, 
and write = The index a takes only two values, 
-|- and — , for positively and negatively charged particles 
respectively. 

The effects of local parity violation cannot be signifi- 
cantly observed in a single event because of the statistical 
fiuctuations in the large number of particles, which are 
not affected by the T'-violating fields. The average of 
over many events, (aa), must be zero. The observation 
of the effect is possible only via correlations, e.g. mea- 
suring [aaap] with the average taken over all events in 
a given event sample. The correlator {aaap) is, however, 
a T'-even quantity, and an experimental measurement of 
this quantity may contain contributions from effects un- 
related to parity violation. The correlator (aaap) can be 
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in principle evaluated via measuring (sin Z\0ct sinZ\(/>^) 
with the average in the last expression taken over all 
pairs of particles of a given type from a kinematic region 
under study and then over all events. The problem is 
that this form of correlator contains also a large contri- 
bution from correlations not related to the reaction plane 
orientation (such correlations are not accounted for by 
Eq. [2) which is a single particle distribution). For exam- 
ple, a pair of particles originating from a single jet will 
typically be emitted with a small azimuthal separation. 
These particle pairs will make a positive contribution to 
{sin A(f)a sinA(j)p), even if the parent jets are emitted 
isotropically relative to the reaction plane. Therefore, 
we separate 



(sinZ\0a sin A(j)p) — (0^0^) -I- Bo 



(3) 



where {aaap) is caused by parity violation and -Bout (de- 
fined by this expression) includes all other correlations 
projected onto the direction perpendicular to the reac- 
tion plane ( "out of plane" ) . The effects contributing to 
Bout may be large and are difficult to estimate reliably. 
For this reason, a different correlator was proposed P^ : 



(cos(0„ + (f)f}- 2^Rp)) = 
(cosZ\(/)q. cos A^ip) — {sin A(f>oc sin 

^^1,^) + Bin] ~ [{O'o.O-p) + Bout 

where, similarly to Eq. |3] Bin is defined via: 

(cOSZ\(?!)q cos A4)p) = {vi^aVl^p) + Bir 



(4) 
(5) 

(6) 



The correlator Eq. [4] represents the difference between 
correlations of the projections of the particle transverse 
momentum unit vectors onto an axis in the reaction plane 
and the correlations of the projections onto an axis that is 
out-of-plane or perpendicular to the reaction plane. The 
key advantage of using Eq. [5] is that it removes all the 
correlations among particles a and /3 that are not related 
to the reaction plane orientation [5D1 ■ 

The contribution given by the term (fi.a'^i,^) can be 
neglected because directed flow averages to zero in a ra- 
pidity region symmetric with respect to mid-rapidity, as 
used in this analysis and the contribution due to directed 
flow fluctuations is very small (see Section VII for a quan- 
titative estimate). Equation [5] then implies that by us- 
ing {cos{(j)a + (jip — 2^Bp)) instead of {sin A(j)a sinZi^/j), 
the background to our measurement of {aaap) is now 
not Bout, but [Bout — Bin], where Bin is the contribu- 
tion of the in-plane correlations which are analogous to 
Bout ■ Only the parts of such correlations that depend on 
azimuthal orientation with respect to the reaction plane 
remain as backgrounds. Studies of the various physics 
contributions to [Bout ~ Bin] are discussed in detail in 
Section Eni 

Based on the current theoretical understanding of the 
Chiral Magnetic Effect one might expect the following 
features of the correlator {uaag): 



• Magnitude: The first estimates [13] predicted a 
signal of the order of |a| ~ Q/N^+, where Q = 
0, ±1, ±2, ... is the net topological charge and A^^+ 
is the positive pion multiplicity in one unit of ra- 
pidity - the expected rapidity scale for correlations 
due to topological domains, see below. More ac- 
curate estimates [TS] including the strength of the 
magnetic field and topological domains production 
rates, were found to be close to the same number. 
It corresponds to values of |a| of the order of 10~^ 
for mid-central collisions, and to 10~^ for the cor- 
relator {aaap). 

• Charge combinations: If the particles, after leav- 
ing the domain experience no medium effects (re- 
interaction with other particles in the system), 
one would expect = — a_. Thus, in the ab- 
sence of medium effects, one expects {aj^.a^) — 
(a_a_) = — (a+a_) > 0. If the process oc- 
curs in a dense medium one needs to account for 
correlation modifications due to particle interac- 
tion with the medium [TS]. The effect of these 
modifications is similar to the modification of the 
jet-like two-particle correlations which experience 
strong suppression of the back-to-back correlations: 
(a+a_|_) = (a_a_) » — (a_|_a_). The effect of 
strong radial flow can further modify this relation 
such that the opposite charge correlations can even 
become positive. 

• Centrality dependence: Under the assumption that 
the average size of the P-violating domain does not 
change with centrality, the correlator should follow 
a 1 /N dependence (typical for any kind of correla- 
tions due to clusters; N is the multiplicity) multi- 
plied by a factor accounting for the variation of the 
magnetic field. The latter is difficult to predict re- 
liably at present, other than that it should be zero 
in perfectly central collisions. Thus at large cen- 
tralities the effect should decrease with centrality 
somewhat faster than 1/N. 

• Rapidity dependence: The correlated particles are 
produced in a domain of the order of 1 fm, and it is 
expected that the correlations should have a width 
in Arj = [rja — '>][}[ of the order unity, as is typical 
for hadronic production from clusters [21] . 

• Transverse momentum dependence: Local parity 
violation is non-perturbative in nature and the 
main contribution to the signal should "come from 
particles which have transverse momentum smaller 
than 1 GeV/c" [15]. The actual limits might be 
affected by the radial flow. 

• Beam species dependence: The effect should be pro- 
portional to the square Z'^ of the nuclear charge, 
but the atomic number A dependence is not well 
understood. One qualitative prediction is that the 
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suppression of the back-to-back correlations should 
be smaller in collisions of lighter nuclei. 

• Collision energy dependence: The effect might be 
stronger at lower energies, as the time integral of 
the magnetic field is larger. At the same time, 
the charge separation effect is expected to depend 
strongly on deconfinement and chiral symmetry 
restoration [TS], and the signal might be greatly 
suppressed or completely absent at an energy below 
that at which a quark-gluon plasma can be formed. 

The main systematic uncertainty in application of the 
correlator Eq.[4]to measurements of anisotropics in parti- 
cle production with respect to the reaction plane, is due 
to processes when particles a and /3 are products of a 
cluster (e.g. resonance, jet) decay, and the cluster itself 
exhibits elliptic flow [22 . Detailed discussion of this 
and other effects which could mimic the effect of local 
strong parity violation in experimental measurements is 
presented in section |VII| 

In this paper, we report our measurements of corre- 
lators shown in Eqs. |4][5] and present systematic studies 
of the background effects that affect the measurements. 
Section |ll] discusses the experimental setup, while sec- 
tion mil discusses the observables and the methods for 
estimating the reaction plane angle and corrections for fi- 
nite reaction plane resolution. Sections IV and [V] present 
the data and a discussion of systematic effects that can 
affect the measurements. Our main results, and how they 
systematically change with system size, centrality, parti- 
cle transverse momentum, and separation in rapidity, are 
presented in Section VI Physics backgrounds that can 



mimic the T'-violating effect are discussed in section [VTTl 



II. EXPERIMENTAL SETUP AND DATA 
TAKING 



The data were collected with the STAR detector at 
Brookhaven National Laboratory during the 2004 and 
2005 runs. Collisions of Au-|-Au and Cu-|-Cu beams 
were recorded at ^/s~nn — 200 and 62 GeV incident 
energies; for a total of four beam-energy combinations. 
Charged particle tracks were reconstructed in a cylindri- 
cal Time Projection Chamber (TPC) [21 US]. The TPC 
is a 4.2 m long barrel with a 2 m radius which was op- 
erated in a solenoidal magnetic field of 0.5 T. The TPC 
detects charged tracks with pseudorapidity |?7| < 1.2 and 
Pt > 100 MeV/c with an absolute efficiency that ranges 
from 80 to 90%. The TPC is nearly azimuthally symmet- 
ric and records tracks at all azimuthal angles; however, 
sector boundaries and other regular detector features are 
responsible for an approximately 10% loss of particles due 
to the finite acceptance of the detector. Track merging, 
and other tracking artifacts that depend on track den- 
sity, can cause an additional 0-10% loss of reconstructed 
tracks; so the overall efficiency is typically 85% per event. 



The TPC's pseudorapidity coverage of an event is sup- 
plemented by two cylindrical and azimuthally symmet- 
ric Forward Time Projection Chambers (FTPC). The 
FTPCs are placed in the forward and backward direction 
relative to the main TPC and cover pseudorapidity inter- 
vals 2.7 < I77I < 3.9 26J. In the most forward direction, 
STAR has two Zero Degree Calorimeter - Shower Maxi- 
mum Detectors (ZDC-SMD) [271 HSj which are sensitive 
to the directed flow of neutrons in the beam rapidity re- 
gions. 

A minimum bias trigger was used during data-taking. 
Events with a primary vertex within 30 cm along the 
beam line from the center of the main TPC were se- 
lected for the analysis. Standard STAR software cuts 
were applied to suppress pile-up and other malformed 
events or tracks. The results presented here are based on 
14. 7M Au+Au and 13. 9M Cu-l-Cu events at the center of 
mass energy of a nucleon pair ^/S7vjv=200 GeV, and 2.4M 
Au-I-Au and 6.3M Cu-I-Cu events at ^5]v^=62 GeV. The 
data were taken with the magnetic field in the Full Field 
(FF), and Reverse Full Field (RFF) configurations, with 
the strength of the magnetic field at 0.5 T. The centrality 
of the collision is determined according to the reference 
multiplicity (refMult), which is the recorded multiplicity 
of charged particles in < 0.5 that satisfy specific track 
quality cuts. 

The correlations are reported in the pseudorapidity re- 
gion I77I < 1.0 covered by the main TPC. For this anal- 
ysis, the tracks in the TPC are required to have pt > 
0.15 GeV/c. For the results integrated over transverse 
momentum we also impose an upper cut of < 2 GeV/c. 
Standard STAR track quality cuts are applied: a mini- 
mum of 15 tracking points are required for a track to be 
considered good. The ratio of the number of hit points 
on a track to the maximum possible given the track ge- 
ometry is required to be greater than 0.52 to avoid the 
effects of track splitting. The data with reverse magnetic 
field were used to assess systematic effects, as the biases 
for positive and negative charged particles interchange. 
The final results reported here are averaged over both 
field polarities. 

We use particle identification via specific energy loss 
[dE/dx) in the volume of the TPC to reject electrons 
as a check that the signal we present is determined by 
hadron production. 



III. METHOD 

In practice the reaction plane angle for a given collision 
is not known. In order to evaluate the correlator defined 
in Eq. |4] one estimates the reaction plane with the so- 
called event plane reconstructed from particle azimuthal 
distributions [25]. For the event plane determination 
one can use particles found in the same detector that is 
used to detect particles a and /3 (in our case STAR's main 
TPC) or different detectors (we have used the STAR FT- 
PCs and the ZDC-SMD). The second order event plane 
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(determined by the second harmonic modulation in par- 
ticle distribution) is sufficient for this study. We make 
use of the large elliptic flow measured at RHIC |3D] to 
determine the event plane from particle distributions in 
the main and Forward TPCs. When using the ZDC-SMD 
for event plane reconstruction, the first-order event plane 
can be determined through the measured directed flow of 
spectator neutrons. 

In the three-particle correlation technique, the explicit 
determination of the event plane is not required; in- 
stead, the role of the event plane is played by the third 
particle that enters the correlator with double the az- 
imuth [351 1221 US • Under the assumption that particle 
c is correlated with particles a and /3 only via common 
correlation to the reaction plane, we have: 

(cos(0a + 0^ - 20c)) = (C0s((/ia + (/'/3 " 2lfflp)) t>2,c, (7) 

where U2,c is the elliptic flow value of the particle c. We 
check this assumption by using particles c from different 
detectors and exhibiting different elliptic flow. We also 
study the effect of using only positive or only negative 
particles to determine the event plane and compare the 
results obtained with different field polarities in our esti- 
mates of the systematic uncertainties. All the correlators 
presented in this paper have been calculated by first av- 
eraging over all particles under study in a given event 
and subsequently averaging the results over all events in 
a given event sample. 

The STAR TPCs have quite uniform azimuthal ac- 
ceptance. Nevertheless, TPC sector boundaries, mal- 
functioning electronics, etc., may introduce biases in the 
analysis, in particular as the acceptance for positive and 
negative particles is different. To avoid these effects we 
use a recentering procedure [29' in which we substitute: 
cos — > cos — (cos (/)) and sin </> — > sin — (sin 0) and 
similarly for the second harmonic. The typical values of 
(cos(/)) and (sin^) for the tracks in the main TPC are 
< 0.003, but for high particles and the most central 
collisions could go as high as 1.5%. In the FTPC region, 
the typical correction is of the order of a few percent. 
The validity of the recentering method can be verified by 
calculating three-particle cumulants [501 [31] • 

(cos(0c< + - 2(?!)c)) = 

+ {upvl) (u„) - 2 {u^) {up) (8) 

where we use notations u = e"^ and v = u* = e~**. 

denotes the real part, and double angle brack- 
ets denote cumulants. In the case of perfect acceptance, 
the cumulant ((cos((/)q + 4>p — 2(f>c))) coincides with the 
correlator (cos(0q + — '^4'c))- As can be seen from 
Eq. [8j to account for the acceptance effect it is sufficient 
to perform a recentering procedure. All results presented 
in this paper have been corrected for acceptance effects, 
where applicable, by this method. The cumulant Eq. [8] 
can be calculated directly by correcting the three-particle 



correlator with the corresponding products of two and 
single particle averages. We have compared the results 
obtained by directly calculating cumulants with the re- 
sults obtained by the recentering method and found them 
to be consistent. Because the detector acceptance varied 
during the period of data taking, we perform the cor- 
responding correction run-by-run, separately for positive 
and negative particles, and for each centrality bin. We 
also account for the acceptance dependence on particle 
pseudorapidity and transverse momentum. We consider 
separately the East (77 < 0) and West {rj > 0) FTPCs. 
We have found that the corrections do not depend sig- 
nificantly on the collision vertex position along the beam 
line. 



IV. DETECTOR EFFECT STUDIES 

Figure [3] shows the three-particle correlator, Eq. [7j as 
a function of reference multiplicity in Au-|-Au collisions 
at y/sNN— 200 GeV for two field polarities before the re- 
centering procedure. All three particles are from main 
TPC region, |?7| < 1.0. Figure |4] shows results for the 
same correlator after correction. The correlator has been 
scaled by the reference multiplicity for clarity at high 
centralities, where the absolute values of the signal are 
small. These figures are intended only to illustrate the 
effect of the recentering; for that reason and also to have 
finer binning in centrality, we plot the correlator directly 
versus reference multiplicity. All other results are pre- 
sented as a function of the fraction of the total interac- 
tion cross section (which is calculated taking into account 
the track, event vertex reconstruction, and trigger inef- 
ficiencies). The acceptance effects are most noticeable 
for central collisions, where the signal is small; there is 
a slight difference in results depending on whether the 
third (c) particle is positive or negative and the difference 
changes sign depending on the polarity of the magnetic 
field. This difference disappears after the acceptance cor- 
rection. Results for particles a and /3 being both positive 
or both negative are consistent within statistical errors, 
and later we combine them as same charge correlations. 
As expected for the case when particles a and /3 are corre- 
lated to the particle c only via common correlation to the 
reaction plane, the results do not depend on the charge 
of the particle c. 

The acceptance effects are larger in the average 
correlation, (cos((/>q — 0^)), than in the correlator 
(cos((/>Q. + (pp — 2'Ffip)), because the latter represents the 
difference in correlations projected onto the reaction 
plane and to the direction normal to the reaction plane. 
As the reaction plane is uniformly distributed in azimuth, 
many of the possible acceptance effects average out to 
zero. 

Figure [5] presents the correlator (cos((/)q — i/)^)) for dif- 
ferent charge combinations from the Au-|-Au 200 GeV 
data obtained with FF and RFF magnetic field settings 
as a function of collision centrality. In this figure and 
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FIG. 3: (Color) (cos(<^q + 4>p — '24'c)) as a function of ref- 
erence multiplicity for different charge combinations, before 
corrections for acceptance effects. In the legend the signs indi- 
cate the charge of particles a, /3, and c. The results shown are 
for Au+Au collisions at 200 GeV obtained in (a) the Reversed 
Full Field, and (b) the Full Field configurations. 



later in the paper, the centrality is quantified by the frac- 
tion of the total interaction cross section, with the cen- 
trality bins corresponding to (ordered from most to least 
central) 0-5%, 5-10%, 10-20%, 70-80% of the most 
central collisions. The points are plotted at the middle of 
the bin, not reflecting possible small biases due to higher 
weight of events with larger multiplicity within the bin. 
Before acceptance corrections are applied, (-I-, -|-) correla- 
tions are slightly different from (— , — ) correlations, with 
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FIG. 4: (Color) Same as Fig.jsjafter correction for acceptance 
effects. 



the difference changing sign in different field orientations. 
After the correction, the results from different field po- 
larities coincide with each other. 

We have performed several additional checks to ensure 
that the signal is not due to detector effects. High ac- 
celerator luminosity leads to significant charge buildup in 
the TPC, which leads to distortions in the recorded track 
positions, affecting the reconstructed momenta. We have 
compared the results obtained from the 2002 RHIC run 
(a low luminosity run), with results from 2004-2005 di- 
vided into high and low luminosity events (selection is 
based on ZDC coincidence rate). All three data samples 
yield the same signal within statistical uncertainties. 
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FIG. 5: (Color) (cos((/)a — 0^)) as a function of centrality 
for different charge combinations and FF and RFF configura- 
tions. The data points corresponding to different charge and 
field configurations are slightly shifted in the horizontal di- 
rection with respect to each other for clarity. The error bars 
are statistical. Also shown are model predictions described in 
Section 



The acceptance of the detector depends weakly on the 
position of the event vertex relative to the center of the 
TPC. We applied the acceptance corrections differen- 
tially according to the event vertex position, and explic- 
itly checked the dependence of the signal on the vertex 
position. No dependence has been found. 

The main TPC consists of two parts which are sepa- 
rated by a central membrane. A particle track will occa- 
sionally cross the central membrane, and be separately 
reconstructed in each half-barrel of the TPC. These two 
track parts can be displaced one with respect to the other. 
In order to check that this effect does not contribute to 
the signal, we calculated the correlator using only tracks 
that do not cross the membrane. Taking into account 
the signal dependence on the track separation in pseudo- 
rapidity, the observed signal was found to be consistent 
with the signal obtained without such a requirement. 

Tracks in the TPC are characterized by the distance 
of closest approach (dca), the distance between the pro- 
jection of the track, and the event vertex. Particles orig- 
inating from weak decays (yl, K^, etc.) can have larger 
dcas than the direct primary particles we are studying. 
We compared the results obtained with a cut dca < 1 cm 
to those of dca < 3 cm, and found only negligible differ- 
ences with a somewhat larger signal (of the order of the 
statistical error) for tracks with dca < 1 cm. 

The correlator used in this analysis is the difference be- 
tween the correlations projected onto the reaction plane 
and the correlations projected onto the direction normal 



to the reaction plane. The correlator calculated by pro- 
jecting onto an axis rotated by 7r/4 relative to the reac- 
tion plane should only be non-zero due to detector ef- 
fects. We have explicitly calculated the correlator in this 
rotated frame and found it to be zero within statistical 
error. 

Figure [6]ja) compares the three-particle correlations 
obtained for different charge combinations, as a function 
of centrality, when the third particle is selected from the 
main TPC with when it is selected from the Forward 
TPCs. Assuming that the second harmonic of the third 
particle is correlated with the first harmonic of the first 
two particles via a common correlation to the reaction 
plane, the correlator should then be proportional to the 
elliptic flow of the third particle. On average, the elliptic 
flow in the FTPC region is significantly smaller than that 
in the TPC region j32| . explaining the different magni- 
tudes of the three-particle correlations shown in Fig. |6ja) . 

Figure [6]Jb) shows the three-particle correlator after it 
has been divided by V2 of the third particle according 
to Eq. [Tj Resulting signals are in very good agreement 
in the two cases. In this and subsequent plots, for the 
elliptic flow of particle c in the main TPC region we use 
estimates obtained from the correlations of particles in 
the main TPC region, \r]\ < 1.0, with particles in the 
FTPC, 2.7 < I77I < 3.9. These estimates are less affected 
by non-flow effects, compared to elliptic flow derived from 
two-particle correlations with both particles taken from 
the main TPC. 

The shaded band in Fig. |6f^b) and the subsequent fig- 
ures illustrate the systematic change in the results that 
occur when different estimates of the elliptic flow are 
used. The upper (in magnitude) limit is obtained with 
flow from four-particle correlations and the lower limit 
from the two-particle cumulant method. All elliptic flow 
data have been taken from Ref. [32j[33] ^. Four-particle 
cumulant values are not available for all collision systems 
and energies studied here. Therefore in Figs. [7]^ we 
plot systematic upper limits obtained with extrapolation 
of available data assuming that the measurements with 
FTPC suppress only 50% of the non-flow contribution. 
The magnitude of the elliptic flow in the FTPC region 
was estimated from correlations between particles in the 
East and West FTPCs. Section |V] has further details on 
the systematic uncertainties associated with different V2 
estimates. 

Results obtained with the event plane reconstructed 
with ZDC-SMD are consistent with those shown in 
Fig. |6|^b), though the statistical errors on ZDC-SMD re- 



in Ref. 1321 1331 an estimate of elliptic flow in the main TPC 
region, \r]\ < 1.0, obtained from correlations of particles in this 
region with those in FTPCs was denoted as ti2{FTPC}; an es- 
timate from two-particle correlations with both particles in the 
main TPC as V2{2}. Elliptic flow from four-particle correlations, 
denoted as );2{4}, is considered to be least affected by non-flow 
effects. For a review of flow measurements, see I18| . 
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FIG. 6: (Color online) (a) A comparison of the correlations 
obtained by selecting the third particle from the main TPC 
or from the Forward TPCs. (b) The results after scaling by 
the flow of the third particle. The shaded areas represent the 
uncertainty from V2,c scaling (see text for details). In both 
panels, the TPC and FTPC points are shifted horizontally 
relative to one another for clarity purposes. The error bars 
are statistical. 



suits are about 5 times larger because the (second order) 
reaction plane resolution from ZDC-SMD is worse. 

Figure |6][b) shows very good agreement between the 
same charge correlations obtained with the third particle 
in the TPC and FTPC regions, which supports for such 
correlations the assumption {cos{(j)a -\- (t)p — 2'Pjip)) w 
(cos((/)q + <^/3 — 2(/)c)) /v2,c- The opposite charge correla- 
tions are small in magnitude and it is difficult to conclude 



on validity of the assumption for such correlations based 
only on results presented in Fig. |6][b). Similarly, in the 
most peripheral collisions, the statistical errors are large, 
which also prohibits making a definite conclusion. 



V. SYSTEMATIC UNCERTAINTIES 

There is one class of uncertainties, related to the ques- 
tion of factorization of Eq. [7j which would arise if the 
events contained a large number of correlated groups of 
particles such as minijets. Even if these "clusters" were 
produced isotropically in azimuth, they might contribute 
to our observable through correlations between the par- 
ticles used to determine the reaction plane (particle c in 
Eq. [t]) and the particles (a, (3) used to measure the sig- 
nal. We consider this effect in detail in Section fVIII As 
will be shown there, in Cu-|-Cu and peripheral Au-|-Au 
collisions this effect could cause opposite charge correla- 
tions of the sign and magnitude we observe, but does not 
produce the same charge correlations. 

We proceed with discussion of the results assuming 

(C0s((/)a +(f>fS~ 'i.l'Rp)) = (C0s((^a + (/"/S " 20c)) /w2,c but 

indicate in all plots the HIJING [37 (default, quenching- 
off settings) three-particle correlation results. The latter 
can be considered as an estimate of the systematic un- 
certainty from correlations not related to the reaction 
plane. In future high statistic measurements such un- 
certainty can be decreased by taking particle c from a 
rapidity region separated from particles a and /3. 

One dominant systematic uncertainty in the correlator 
(cos((/)q + 0^ — 2'Fjip)) is due to uncertainty in the ellip- 
tic ffow measurements of the particle used to determine 
the reaction plane. This contributes a fractional uncer- 
tainty, on average of the order of 15% and somewhat 
larger in most peripheral and most central collisions [32j . 

From comparison of the results obtained in different 
field configurations, and other studies presented in Sec- 



tion IV we conclude that after acceptance corrections 
are performed, the remaining systematic uncertainties 
in three-particle correlations due to detector effects are 
comparable to or smaller than the statistical errors. 

We have performed an additional study to estimate 
the size of possible error caused by acceptance effects be- 
fore and after the recentering correction is applied: we 
have run simulations in which tracks were generated us- 
ing realistic single particle distributions but having no 
correlation except due to elliptic flow. An efficiency loss 
is introduced similar to that of the STAR detector as a 
function of azimuth, transverse momentum, and particle 
charge. We then study the effect of distorting the effi- 
ciency in additional and more extreme ways. In all of 
these cases, after the recentering correction is applied, 
the value of (cos((/>q, + (/)^ — 20c)) is zero for all centrali- 
ties within the statistical precision of the study, which is 
about 3 X 10~^ for the most peripheral bin and decreases 
to less than 10~^ for the most central bin. This is many 
times smaller than the measured signal for all centralities 
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in all cases. 

Errors in measuring the magnitude of particle mo- 
menta make negligible contributions to the correlator 
used in this analysis, which uses only measured azimuthal 
angles. It is therefore robust against many systematic er- 
rors which are commonly encountered in the analyses of 
the STAR data (space charge distortion errors leading to 
momentum biases, etc.). 

Theoretical treatments of the correlator defined in 
Eq. [4]were developed with charged hadrons in mind. By 
using cuts (based on specific energy loss) to suppress the 
presence of electrons in our sample, we have verified that 
this bias is also smaller than the statistical errors. 



VI. RESULTS 
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Final results presented in this section have been ob- 
tained with three-particle correlation using Eq.[7|with aU 
three particles from the pseudorapidity region |?7| < 1.0. 
Figure [7| presents the correlator {cos{(f>a + 0/3 — 2'Fiip)) 
for Au-I-Au and Cu-|-Cu collisions at ^/sn~n— 200 GeV. 
Positive-positive and negative-negative correlations are 
found to be the same within statistical errors, see 
Fig. Qb), and are combined together as same-charge 
correlations. Opposite-charge correlations are relatively 
smaller than same-charge correlations, in agreement with 
possible suppression of the back-to-back correlations dis- 
cussed in the introduction. The correlations in Cu+Cu 
collisions, shown by open symbols, appear to be larger 
than the correlations in Au+Au for the same centrality 
of the collision. One reason for this difference may be the 
difference in number of participants (or charge multiplic- 
ity) in Au-I-Au and Cu-|-Cu collisions at the same cen- 
trality. The signal is expected to have a 1/7V dependence, 
and at the same centrality of the collision the multiplicity 
is smaller in Cu+Cu collisions than in Au+Au. The dif- 
ference in magnitude between same and opposite charge 
correlations is considerably smaller in Cu+Cu than in 
Au+Au, qualitatively in agreement with the scenario of 
stronger suppression of the back-to-back correlations in 
Au+Au collisions. In Fig. [7] and below, error bars indi- 
cate statistical uncertainties. The shaded bands show the 
systematic uncertainty associated with measurements of 
elliptic flow which are used to rescale the three-particle 
correlator. In this section we assume the factorization 
of correlator Eq. [7] The possible error due to this as- 
sumption - which may be large for peripheral bins in 
the opposite charge correlation - is denoted by the thick 
lines in F ig. [7| and subsequent figures and is explained 
in Section |Vnj Other systematic uncertainties have been 
discussed in section FVl 

Figure [8] shows results for collisions at 
y/sNN=62A GeV. The signal is similar in magni- 
tude, with slightly more pronounced opposite-charge 
correlations compared to those at ^/snn=200 GeV. This 
is consistent with weaker suppression of opposite-charge 
correlations in the less dense 62 GeV system. 



FIG. 7: (Color online) {cos{(j)a + <Pi3 ~ 2^i{p)) in Au+Au and 
Cu+Cu collisions at ^/snn = 200 GeV calculated using Eq.[7| 
The error-bars show the statistical errors. The shaded area 
reflects the uncertainty in the elliptic flow values used in cal- 
culations, with lower (in magnitude) limit obtained with ellip- 
tic flow from two-particle correlations and upper limit from 
four-particle cumulants. For details, see Section |IV[ Thick 
solid (Au+Au) and dashed (Cu+Cu) lines represent possi- 
ble non-reaction-plane dependent contribution from many- 
particle clusters as estimated by HIJING, see Section [VII A| 
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FIG. 9; (Color online) Au+Au and Cu+Cu collisions at 
V'SiVJV=200 GeV. The correlations are scaled with the num- 
ber of participants and are plotted as function of (a) centrality 
and (b) number of participants. The error-bars indicate the 
statistical errors. The shaded area reflects the uncertainty 
in the elliptic flow values used in calculations. For details, 



see Section IV Thick solid (Au-fAu) and dashed (Cu-f Cu) 
lines represent possible non-reaction-plane dependent contri- 
bution from many-particle clusters as estimated by HIJING, 
see Section rVilAI 



The correlations are weaker in more central collisions 
compared to more peripheral collisions, which partially 
can be attributed to dilution of correlations which occurs 
in the case of particle production from multiple sources. 
To compensate for this effect and to present a more com- 
plete picture of the centrality dependence, we show in 
Fig. [9] results multiplied by the number of participants. 



The number of nucleon participants is estimated from a 
Monte-Carlo Glauber model |35j. Figure 9]ja) presents 
the results as a function of centrality, and Fig. [9]Jb) as 



function of N, 



part • 



Smaller correlations in most cen- 



tral collisions are expected in the parity violation picture 
as the magnetic field weakens. The same and opposite 
charge correlations clearly exhibit very different behav- 
ior. Figure [9f^a) demonstrates that the same-charge cor- 
relations show similar centrality dependencies, as would 
be expected if the geometry of the collision is important. 
The opposite-charge correlations in Au-|-Au and Cu-|-Cu 
collisions are found to be close at similar values of Npart 
in rough qualitative agreement with the picture in which 
their values are mostly determined by the suppression of 
back-to-back correlations. 

Figure [TO] shows the dependence of the signal on the 
difference in pseudorapidities of two particles, Arj = 
\Va — Vi3\j for 30-50% and 10-30% centralities. The signal 
has a typical hadronic width of about one unit of pseudo- 
rapidity. The dependence on |?7q,— ?7^| has been calculated 
for all charged tracks with 0.15 < pt < 2.0 GeV/c. Fig- 
ure 11 shows the dependence of the signal on the sum of 
the transverse momentum (magnitudes) of the two parti- 
cles for these same centralities. Results presented is this 
figure have no upper pt cut. We do not observe the sig- 
nal concentration in the low pt region as naively might 
be expected for T'-violation effects. 

Figure [12] shows the dependence of the signal on the 
difference in the magnitudes of the two particle trans- 
verse momenta. We find that the correlation depends 
very weakly on \pt,a — Pt,i3\- This excludes quantum in- 
terference (HBT) or Coulomb effects as possible expla- 
nations for the signal. There are no specific theoretical 
predictions on this dependence for the chiral magnetic 
effect, though naively one expects that the signal should 
not extend to large values of \pt,a — Pt,i3\- 

Finally, the ZDC-SMD detector has good first-order 
reaction-plane resolution. For mid-central collisions the 
resolution, (cos('f'i — ^rp)), is of the order of 0.35-0.4. 
The ZDC-SMD allows us to test the first-order [V-oAA) 
effect of the charge separation along the system orbital 
momentum, which would correspond to (0^) ^ 0. In 
theory this is possible only if the vacuum 9 ^ Q. Our 
measurements are consistent with zero, averaged over all 
centralities (a+) = (-0.1± 1.0) • 10"* and (a_) = (-1.0± 
1.0) • 10-4. 



VII. PHYSICS BACKGROUNDS 

A. Reaction-Plane Independent Background 

Reaction-plane independent background is caused by 
three (or more) particle clusters which affect the factor- 
ization of Eq. [7] With future high statistics data sets, it 
will be possible to reduce such backgrounds significantly 
by determining the reaction plane using particles far re- 
mote in rapidity from the signal particles. 
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FIG. 10: (Color online) Au+Au at 200 GeV. The correlations 
dependence on pseudorapidity separation Arj = — 77/3I for 
(a) centrality 30-50%, and (b) centrality 10-30%. The shaded 
band indicates uncertainty associated with V2 measurements 
and has been calculated using two- and four-particle cumulant 
results as the limits. 



In order to estimate possible contribution to the three- 
particle correlator of effects not related to the reaction 
plane orientation we use the HIJING [Mj event generator 
which is based on the minijet picture of heavy-ion colli- 
sions. For all HIJING results presented in this paper we 
use default, quenching-off setting. Figure [13] presents the 
results for the three-particle correlator, (cos(^a + 0^9 — 
2(f>c)), measured in Au-|-Au and Cu-|-Gu coUisions as a 
function of centrality together with HIJING results for 
the correlations among three particles from many parti- 
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FIG. 11: (Color online) Au-|-Au at 200 GeV. The correla- 
tions dependence on {pt,a +pt,p)/'2 for (a) centrality 30-50%, 
and (b) centrality 10-30%. The shaded band has the same 
meaning as that in Fig. |10| 



cle clusters. In this figure the most central points cor- 
respond to centrality 0-5% and the most peripheral to 
60-70% for Au-HAu coUisions and 50-60% for Cu-hCu 
collisions. The correlator is scaled with number of par- 
ticipants for clarity at large centralities, where the signal 
is small in magnitude. The correlations are shown as a 
function of the number of participants because this gives 
very similar HIJING results for Au-f Au and Cu-|-Cu col- 
lisions, implying a dependence only on the charged par- 
ticle rapidity density. We have separately checked that 
HIJING results scale as A^~^ as expected for contribu- 



tions from many particle clusters. Figure 13 shows that 



14 



,x10"- 



0. 

1.5 



+ « 0.5 



§ 

^ -0.5 



-1 
-1.5 
-2 



STAR AuAu 200 GeV 

Centrality 30-50% 
— • — same charge 
— ■ — opp charge 



, *Ar^^^^^^^^^^^^^^^^^^^^^ i^^^'f^ 



_l I I I L_ 



_l I I I L_ 



0.5 1 1.5 2 

(a) 



,x10"- 



0. 

1.5 

IF 

-e- 

+ a 0.5 r 



-1 r 
-1.5 r 

-2p- 



T" 



STAR AuAu 200 GeV 

Centrality 10-30% 
— • — same charge 
— ■ — opp charge 



_l I I I L_ 



_l I I I L_ 



0.5 1 1.5 2 

IPt,a-Pt,,l (G«V/C) 

(b) 

FIG. 12: (Color online) Au+Au at 200 GeV. The correlations 
dependence on \pt,a — Pt,/?! for (a) centrality 30-50%, and (b) 
centrality 10-30%. The shaded band has the same meaning 
as that in Fig. [lO] 



if this minijet picture in HIJING is correct, in peripheral 
collisions the entire opposite charge signal may be dom- 
inated by contributions from clusters not related to the 
reaction plane orientation. The same charge correlations 
in HIJING are significantly smaller in magnitude than in 
data and have opposite sign. HIJING results for three- 
particle correlations among three particles all of the same 
charge are consistent with zero in sharp contrast to the 
data shown in Fig. |4] 

We have also studied such reaction-plane independent 
backgrounds using the event generator UrQMD and 
found that the predicted contributions to both opposite- 



charge and same-charge correlations are at least a factor 
of two lower than those predicted by HIJING. 
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FIG. 13: (Color online) Three-particle correlator in Au+Au 
and Cu-|-Cu collisions compared to HIJING calculations 
shown as thick lines. All three particles are taken in the main 
TPC region, jjyj < 1.0. The correlator has been scaled with 
number of participants and are plotted versus number of par- 
ticipants. In this representation HIJING results for Au-|-Au 
and Cu+Cu collisions coincide in the region of overlap. 



B. Reaction-Plane Dependent Background 



Unlike those discussed in Sec. |VllXl reaction-plane de- 
pendent physics backgrounds can not be suppressed by 
better methods of determining the reaction plane. 

The correlator (cos(0q + 0/3 — "^^rp)) is a P-even ob- 
servable and can exhibit a non-zero signal for effects not 
related to T'-violation. Among those are processes in 
which particles a and /3 are products of a cluster (e.g. 
resonance, jet, di-jets) decay, and the cluster itself ex- 
hibits elliptic flow [321 [22 or decays (fragments) differ- 
ently when emitted in-plane compared to out-of-plane. 

If "flowing clusters" are the only contribution to the 
correlator, we can write: 



(COS(0, 
Adust (cos((0Q -f 0^ 
= Adust (cos(0Q 



'2'4>dust) + 2(0ciMst ~ ^Rp)))dust 
\- (j)/} — 2(j3dust))dust V2,dust, (9) 



where {■■■)dust indicates that the average is performed 
only over pairs consisting of two daughters from the same 
cluster and the resulting normalization factor is Adust = 
N au^N pairs /N pairs . Equatiou [ol assumcs that there is 

eijenf dust euent J_J 

no reaction plane dependence oicos(0q + 0/3 ^ '^4'dust)- 
The term (cos(0a + 0/3 — '24'dust)) is a measure of the 
azimuthal correlations of decay products with respect to 
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the cluster azimuth, while V2,ciust is cluster elliptic flow. 
In the case of resonance decays, (cos(0a + 0^ — 2(f)res)) 
is zero if the resonance is at rest, and becomes non-zero 
only due to resonance motion. Estimates of the con- 
tribution of "flowing resonances" , based on Eq. [9] and 
reasonable values of resonance abundances and values of 
elliptic flow, indicate that they should not produce a fake 
signal. Given the relative scarcity of parents decaying to 
two same-charge daughters, a much smaller magnitude 
is expected for same charge than opposite charge corre- 
lations from this source. Kinematic studies demonstrate 
that it is very difficult for the sign of the correlations 
observed in the data to be created in the same-charge 
correlations without postulating a negative value of V2 
for the resonances or particles from cluster decays. 

To study the contribution from resonances in greater 
detail we have carried out simulations using the MEVSIM 
event generator [36]. MEVSIM generates particles ac- 
cording to the single particle momentum distributions 
measured at RHIC. The only correlations included are 
an overall bulk elliptic flow pattern and correlations be- 
tween daughters of the same resonance decay (resonances 
included are 0, Z\, p, w, and K*). MEVSIM simulation 



results are shown as solid squares in Fig. 14 the oppo- 
site charge correlations are larger than what is seen in the 
data, while the same charge correlations are far smaller 
in magnitude and of the wrong sign to match parity vio- 
lation correlations. We conclude that resonances are not 
responsible for the observed signal. 

In addition to contributing to reaction-plane indepen- 
dent background as discussed in Section |VII A[ jets 
are another potential source of reaction-plane dependent 
background since their properties may vary with respect 
to the reaction plane. For those jets in a heavy-ion event 
which include a charged particle of sufficiently high pt to 
act as a trigger particle for a jet analysis, we may esti- 
mate the contribution using the results of previous STAR 
jet studies f37H39j . With trigger transverse momentum 
values that allow such analysis {pt > 3 GeV/c) the contri- 
bution to (cos(0c( + 4>i3 — is roughly two orders of 
magnitude below the same charge signal shown in Fig. [7] 
To extend the study of jet contributions to lower momen- 
tum, we rely on event generators (in particular, HIJING) 
calculations. 

Several correlation measurements from RHIC 1301 SI] 
and earlier measurements at ISR (see review |2T) ) indi- 
cate that cluster formation plays an important role in 
multiparticle production at high energies. These clus- 
ters, with a size inferred in KJ] to be 2.5-3 charged 
particles per cluster, may account for production of a sig- 
nificant fraction of all particles. Because we have limited 
information about the nature of these clusters we do not 
make an estimate of their contribution to the observed 
correlations. Our studies indicate that in order to fake 
the same-charge correlations observed in the data, there 
should be several types of clusters with some of them 
having negative values of elliptic flow. It is hoped that 
with a better understanding of the cause and properties 



(including charge dependence and V2) of such clusters, a 
clearer statement can be made regarding their contribu- 
tions. 

We have also run simulations with several j>+p and 
Au-I-Au event generators. With PYTHIA [42, p-l-p 
events we find that the correlations in (cos(<^q, — (pp)) 
are significantly smaller than those seen in Au-|-Au data 
when scaled by 1/A^, and are similar for all charge com- 
binations. We add modulation with respect to the reac- 
tion plane by adding V2 through angular correlations or 
strong (elliptically modulated) radial flow. This way we 
create non-zero values for (cos(0a + (pp — 2^j^p)), albeit 
with correlations different from the data, being always 
positive and similar in magnitude for all charge combi- 
nations. 



Figure 14 shows results for (reaction-plane depen- 
dent) physics backgrounds to (cos((/>a -\- — 2!Z'/jp)) cal- 
culated with 200 GeV Au+Au events from the event gen- 
erators UrQMD 43J and HIJING [34 . Because the mod- 
ulation of dN/d(j) with respect to the reaction plane is 
smaller in HIJING than seen in RHIC data, we also run 
HIJING with an added "afterburner" which adds ellip- 
tic flow using as input V2 values consistent with STAR 
measurements at the given centrality. Elliptic flow is in- 
troduced by the "shifting" method [55] j which preserves 
other correlations that exist in the model. Figure [M] 
shows that no generator gives qualitative agreement with 
the data; the model values of (cos(<^a + <l)j3 — 2^fip)) are 
signiflcantly smaller in magnitude than what is seen in 
the data, and the correlations calculated in these mod- 
els tend to be very similar for same and opposite charge 
correlations. 

These models do not match the correlations for 
(cos((/)q, — (pp)) that are seen in the data either, as shown 
in Fig. [5] HIJING predicts very similar same and 
opposite charge correlations that are much smaller in 
magnitude than seen in the data. UrQMD overesti- 
mates the same charge correlations. It predicts oppo- 
site charge correlations that are much smaller in mag- 
nitude and opposite in sign from the data. This points 
to the need for better modeling of two-particle correla- 
tions to give quantitatively meaningful comparisons for 

(cOs((/)„ -h</)/3 -2!fijp)). 

In Fig. [T4| we connect UrQMD points by dashed lines 
to illustrate that the "reference line" for strong parity 
correlations might be not at zero. In this particular case 
of UrQMD, both same and opposite-charge correlations 
have values below zero. Note that the same-charge cor- 
relations sit somewhat above the opposite-charge corre- 
lations, opposite to the expectation from local parity vi- 
olation. 

Directed flow, which on average is zero in a symmetric 
pseudorapidity interval, can contribute to the correlator 
(cos((/)q, + (pp — 2'Ffip)) via flow fluctuations. This effect 
is of the opposite sign, see Eq. [5| and is similar for differ- 
ent charge combinations unlike the signal. If one assumes 
that the amplitude of the fluctuations is of the same or- 
der of magnitude as the maximum directed flow in the 
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FIG. 14: (Color) {cos((^q + <^,9 - 2<Prp)) calculated for 
200 GeV Au+Au events with event generators HIJING (with 
and without an "elliptic flow afterburner"), UrQMD, and 
MEVSIM. Blue symbols mark opposite-charge correlations, 
and red are same-charge. Solid stars represent the values 
from the data to facilitate comparison. Acceptance cuts of 
0.15 < pt < 2 GeV/c and \ri\ < 1.0 were used in all cases. 
For MEVSIM, HIJING, and UrQMD points the true reaction 
plane from the generated event was used for •Z'bp. Thick 
solid lighter colored lines represent possible non-reaction- 
plane dependent contribution from many-particle clusters as 
estimated by HIJING and discussed in Section |VII A[ Cor- 
responding estimates from UrQMD are about factor of two 
smaller. 



collisions has been presented for Au-|-Au and Cu-|-Cu 
collisions at y's]v]v=200 and 62 GeV. The results are 
reported for different particle charge combinations as 
a function of collision centrality, particle separation 
in pseudorapidity, and particle transverse momentum. 
Qualitatively the results agree with the magnitude and 
gross features of the theoretical predictions for local V- 
violation in heavy-ion collisions, except that the signal 
persists to higher transverse momenta than expected [H] . 
The particular observable used in our analysis is P-even 
and might be sensitive to non-parity-violating effects. So 
far, with the systematics checks discussed in this pa- 
per, we have not identified effects that would explain 
the observed same-charge correlations. The observed 
signal cannot be described by the background models 
that we have studied (HIJING, HIJING-I-U2, UrQMD, 
MEVSIM), which span a broad range of hadronic physics. 

A number of future experiments and analyses are nat- 
urally suggested by these results. One of them is the 
study of the correlation dependence on the energy of the 
colliding ions. The charge separation effect is expected 
to depend strongly on the formation of a quark-gluon 
plasma [TS], and the signal might be greatly suppressed 
or completely absent at an energy below that at which a 
quark-gluon plasma can be formed. 

Improved theoretical calculations of the expected sig- 
nal and potential physics backgrounds in high energy 
heavy ion collisions are essential to understand whether 
or not the observed signal is due to local strong parity 
violation, and to further experimental study of this phe- 
nomenon. 



pseudorapidity interval under study, then the flow fluc- 
tuation contribution is no more than 10~^ for centrality 
30-60%, significantly smaller than the observed signal. 

Global polarization of hyperons [IH |45], the phe- 
nomenon of the polarization of secondary produced par- 
ticles along the direction of the system's angular momen- 
tum, may also contribute to the correlator Eq. |4] via V- 
odd weak decays. This effect could lead to a charge asym- 
metry with respect to the reaction plane, which is always 
pointing in the same direction relative to the orientation 
of the angular momentum. Our main analysis based on 
the reaction plane reconstructed from the elliptic flow 
does not distinguish the direction of the angular momen- 
tum, and is susceptible to this effect. But as we pointed 



out in section VI our measurement of charge separa- 
tion along the system orbital angular momentum is zero 
based on the flrst-order reaction plane reconstructed in 
the ZDC-SMD. Global polarization has also been found 
to be consistent with zero. Pa. a < 0-02 [46] . 

VIII. SUMMARY 
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